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DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
COMMON, as to the accuracy and functioning of the program 
and related program material, nor shall the fact of distribution 
constitute any such warranty, and no responsibility is assumed 
by the contributor or COMMON, in connection therewith. 



COMMON USERS GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter, ink or pencil) 

Program No._ Date 

Program Name: _ 

1. Does the abstract adequately describe what the program is and what Yes TSTo_ 

it does ? 

Comment 

2. Does the program do what the abstract says ? Yes No_ 

C om m ent . _______ 

3. Is the description clear, understandable, and adequate? Yes No_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No_ 

C om m ent 

Are the Sense Switch options adequately described (if applicable) ? Yes No_ 
Are the mnemonic labels identified or sufficiently understandable? Yes No_ 

Comment 



5. Does the source program compile satisfactorily (if applicable)? Yes_ No_ 

Comment 

6. Does the object program run satisfactorily? Yes_ No_ 

Comment 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

Comment 

8. Does the Program meet the minimal standards of COMMON? Yes No_ 

Comment 



9. Were all necessary parts of the program received? Yes 

Comment _ 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 

Company 

Mr. Richard L. Pratt . 

Data Corporation Address 

7500 Old Xenia Pike 



No- 



Dayton, Ohio 45432 Users Group Code ' 

THIS REVIEW FORM _3 PART OF THE COMMON ORGANIZATION'S PROGRAM REVIEW AND 
EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVITED TO PARTICIPATE 
IN THIS EVA LUATION. 
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MOSES 

MAXIMUM ONE-PASS SIMULTANEOUS 
EQUATIONS SOLUTION 



RONALD HOLLMEIER 
PIONEER SERVICE & ENGINEERING CO. 
2 NORTH RIVERSIDE PIAZA 
CHICAGO, ILLINOIS 60606 



3111 



FEBRUARY 11, i960 



Modifications or revisions to this program, as they occur, will be announced 
in the appropriate Catalog of Programs for IBM Data Processing Systems. When 
such an announcement occurs, users should order a complete new program from 
the Program Information Department. 
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1620 USERS GROUP UBRARY 
PROGRAM ABSTRACT 



1. TITLE (If subroutine, state in Title): MOSES - Maximum One-Bass Simultaneous 
Equations Solution ' ■ Subject Classification: ^.o 

2. Author; Organization: Ronald Hollmeier. Pioneer Service & Engineering Co. 

Date: February 11, 1966 Users Group Membership Code: 3111 

3. Direct Inquiries to Name: Ronald Hollmeier, Pioneer Service & Engineering Co.. 

Chicago ? Illinois ; Phone: 312-CE6-7272 

4. Description/Purpose: (5. Method; 6. Restriction/Range; When Applicable): 

Solves sets of linear simultaneous equations by gaussian elimination. Pivoting for 
size to promote accuracy and eliminate zero diagonal terms Is automatic. Equation 
sets of rank 2 thru hi can be solved. One constant vector per equation set is allowed . 



7. Specifications (Check or fill in appropriate spaces): 

a. Storage used by program: 20 k core 

b. Equipment required by program: Car d x ; Magnetic Tape ; Number of 

Drives ; Paper Tape ; Disk File_ ; Number of Drives ; TNS, 

TNF, MF ; Auto divide ; Indirect addressing ; Floating Point 

Hardware ; 1620 Model I x ; Model II ; 1443 Printer ; Index 

Registers ; Binary Capabilities ; Other (specify) 

Can program be used, on lesser machine? no Specify which requirements 
can be easily removed^ . 

c. Programmed in: Fortran without Forma t x ; Fortran with Format ; 

Fortran BE ; Other Fortran (specify) ; SPS (specify assembler used) 

; Other (specify) 

d. Type of Program: Mainline, complete ; Subroutine ; If subroutine, for 

use with SPS (specify type of SPS) ; Fortran (specify type of Fortran 

; Other (specify) ; . 

8. Additional Remarks: 
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DECK LABELLING KEY 



DECK NO. 1 

Source deck 51 cards 

Sequence No. 301^0020 thru 301^0520 



DECK NO. 2 

Object deck 389 cards 
Sequence No. 0001 thru 0389 



DECK NO. 3 

Sample problem input, data if 3 cards 
Sequence No. 001 thru 0^3 



DECK NO. h 

Sample problem results 6 cards 
Sequence No. 001 thru 006 




MOSES 

MAXIMUM ONE-PASS SIMULTANEOUS 
EQUATIONS SOLUTION 



FEBRUARY 11, I966 



RONALD HOLLMEIER 
PIONEER SERVICE & ENGINEERING CO. 
2 NORTH RIVERSIDE PLAZA 
CHICAGO, ILLINOIS 60606 
312-CE6-7272 



3111 



Moses was written to provide a fast and flexible one pass program for solving 
simultaneous equations sets. The program was written to use the Fortran with- 
out Format compiler because it provides a fast object program with sufficient 
input-output subroutines for ease of operation. An "ABS" subroutine is used 
and therefore must be among the available subroutines if the program is modified 
and recompiled. This is mentioned because the "ABS" routine was not part of our 
original Fortran without Format subroutine set. The object program and subroutines 
occupy core locations 00000 thru 10026. The symbol table is in lOHO thru 19999 

The program will solve any nonsingular set of N equations with N unknown having 
from 2 thru kl equations. Using simple Gaussian elimination each equation is 
completely processed after it is read-in. Pivoting is on the largest term in 
each row. This is done to maintain as much accuracy as possible with the eight 
digit limit on mantissa length, and also eliminates any zero diagonal terms. 
Only terms necessary to produce the solution after the elimination process is 
completed are saved. 

Solutions computed using this program have been checked against those produced 
by Simultaneous Equations Ala-King program No. 5-0.008. . This was done for 
matrices of rank 6, 2k, 36, and 1^1. Results compared favorably in all cases. 
However, no attempt at error analysis has been made, and the eight digit mantissa 
may limit the usefulness of this program where extreme accuracy is required. 
Program running times, exclusive of loading, range from approximately 25 seconds, 
to 26.25 minutes for 6 and hi equation matrices respectively. 
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Input requirements, in order, are as follows: 

X) N - Number of equations 

2) B(l) — B(N) - Equation constants 

3) A(l) --- A(N*N) - Equation coefficients 

Coefficients are entered row by row. Input is in free form subject only to 
the usual requirements of the Fortran without Format subroutines. 

Output consists of equation solutions, one per card in order X(l) --- X(N) . 

No sense switches are used. The program will work with any 1620 Model I having 
card reader and punch. 
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"SAMPLE PROBLEM INPUT DATA 



-24719.9340000000 
28319.4220000000 
-5407.6971000000 
.0000000000 
.0000000000 
.0000000000 
294.3113200000 
39.4664600000 
-30.9567290000 
2.4579793000 
5.3140150000 
-7.9745280000 
39.46646UUUUU 
373.93868UOOOU 
70.5523170000 
-.1489070000 
-7.0165537000 
-70.9269090000 
-30.9567320000 
70.5523160000 
397.5333700000 
7.7321770000 

81.2899830000 
4.5585743000 
2.4579793000 
-.1489070000 
7.7321770000 

33.3270330000 
-.7021131700 
1.1670051000 
5.3140150000 

-7.0165536000 

81.2899830000 
-.7021131600 

33.2608180000 
2.6584285000 

-7.9745280000 
-70.9269090000 
4.5585744000 
1.1670051000 
2.6584285000 

29.7637160000 



SAMPLE PROBLEM RESULTS 

-121.88338 

208.05722 
-143.17812 

3.5.348009 

378.07925 

449.91809 



PROGRAM LISTING 





DIMENSION K141). A<862>» 


B 1 4 1 ) 


— 






READ* N 








■?.y A4vV-2.y. 




DO 20 I = 1* 


N 






JU l*»UU*rU 


20 


READ* B(I) 









30140050 




LL = 








3Ul*H)GoU 




NN = N - 1 








A.4y.v f..v„ 




DO 140 M = 1 » 


N 






30 140080 




L = LL + 1 








Ort 1 A Art on 




LL = LL + N 












DO 40 J = Li 


LL 




- 






READ* A(J) 








3 rt 1 A rt i tr\ 




IFIL - 1) 90. 


90. 50 






30140130 


50 


JJ = 








30140140 




MN = M - 1 








3014015J5 




DO 80 I = h 


MN 






30140160 




J = K ( I ) + L 












DUM = A(J) 












A(J) = AIL) 








im A rt 1 On 
301**0i70 




LL = LL - 1 








30 1402 00 




DO 70 J = Li 


LL 










JJ = JJ + 1 










70 


A ( J ) = A( J + l ) 


- DUM*A(JJ) 






art i Art 9 


80 


B(M) = B(M) - 


DUM*B( I ) 






0. n 1 Art O Art 


90 


DUM = 0. 








■art i /.rt 9 c>Ci 




DO 110 J = L » 


LL 






•art 1 Art 




IFIDUM - ABS ( A ( J ) ) ) 100. 








1 


DUM = ABS ( A ( J ) ) 






artlArtOftn 




JJ = J 








a n 1 A rt OQA 




CONT INUE 








30 140300 




DUM = A(JJ) 












A(JJ) = A ( L ) 








30140320 




LL = LL - 1 








30140330 




DO 120 J = Li 


LL 






30140340 


120 


A(J) = A (J+l) /DUM 






3014Q350 




B(M) = B(M)/DUM 






30140360 


140 


KIM) = JJ - L 








30140370 




DO 170 M = 1. 


NN 






30140380 




I = N - M 








30140390 




LL = L - 1 








30140400 




L = L - M 








30140410. 




DUM = B I I ) 








30140420 




J J = K I I ) + I 








30140430 




BII) = B ( J J ) 








30140440 




DO 160 J = L » 


LL 






30140450 




1 = 1 + 1 








30140460 


160 


DUM = DUM - A I J ) *B ( I ) 








170 


BIJJ) = DUM 








30140480 




DO 180 I = 1. 


N 






30140490 


180 


PUNCH* BII) 








30140500 




STOP 








3©14Q_5_1Q 




END 








30140520 




